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Spin-dependent transport measurements through a double quantum dot are a valuable tool for 
detecting both the coherent evolution of the spin state of a single electron as well as the hybridization 
of two-electron spin states. In this paper, we discuss a model that describes the transport cycle in 
this regime, including the effects of an oscillating magnetic field (causing electron spin resonance) 
and the effective nuclear fields on the spin states in the two dots. We numerically calculate the 
current flow due to the induced spin flips via electron spin resonance and we study the detector 
efficiency for a range of parameters. The experimental data are compared with the model and we 
find a reasonable agreement. 



A. Introduction 

Recently, coherent spin rotations of a single electron 
were demonstrated in a double quantum dot device [1]. 
In this system, spin-flips of an electron in the dot were in- 
duced via an oscillating magnetic field (electron spin res- 
onance or ESR) and detected through a spin-dependent 
transition of the electron to another dot, which already 
contained one additional electron. This detection scheme 
is an extension of the proposal for ESR detection in a sin- 
gle quantum dot by Engel and Loss [2] . Briefly, the device 
can be operated (in a spin blockade regime [3] ) such that 
the electron in the left dot can only move to the right dot 
if a spin flip in one of the two dots is induced via ESR. 
From the right dot, the electron exits to the right reser- 
voir and another electron enters the left dot from the left 
reservoir. A continuous repetition of this transition will 
result in a net current flow. 

Compared to the single dot detection scheme [2], us- 
ing the double-dot as the detector has two major advan- 
tages. First, the experiment can be performed at a lower 
static magnetic field and consequently with lower, tech- 
nically less demanding, excitation frequencies. Second, 
the spin detection is rather insensitive to unwanted oscil- 
lating electric fields, because the relevant dot levels can 
be positioned far from the Fermi energies of the leads. 
These electric fields are unavoidably generated together 
with the oscillating magnetic field as well. 

The drawback of the double-dot detector is that spin 
detection is based on the projection in the two-electron 
singlet-triplet basis, while the aim is to detect single spin 
rotations. However, this detection is still possible be- 
cause the electrons in the two dots experience different 
effective nuclear fields. This is due to the hyperfine inter- 
action of the electron spins with the (roughly 10 6 ) nuclear 
spins in the host semiconductor material of each quan- 
tum dot [4-11]. In order to provide more insight in this 
double-dot ESR detection scheme for single spin rota- 
tions, it is necessary to analyze the coherent evolution of 
the two-electron spin states together with the transitions 
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in the transport cycle. 

In this paper, we discuss a model that describes the 
transport cycle in the spin blockade regime while includ- 
ing the coherent coupling between the two dots, and the 
influence of the static and oscillating magnetic field to- 
gether with the effective nuclear fields on the electron 
spin states. The aim is to understand how effectively 
single spin resonance will affect the measured quantity 
in the experiment, namely the current flow in the spin 
blockade regime. The organization of this paper is as 
follows. First, we will explain the transport cycle and 
the mechanism that causes spin blockade. Next, we will 
briefly discuss the static system Hamiltonian and the 
mixing of the two-electron spin states by the effective 
nuclear field. Then we add an oscillating magnetic field 
to this Hamiltonian, that forms -together with the double 
dot tunnelling processes- the basis of the rate equations 
that describe how the density matrix of the two-electron 
spin states evolves in time. The current flow, calculated 
from the steady state solution of the density operator, 
is then analyzed for different coherent coupling values, 
magnitudes of the oscillating magnetic field, in combina- 
tion with different effective nuclear fields in the two dots. 
This provides further insight in the optimal conditions 
for spin-flip detection with a double quantum dot. 



B. Spin blockade 

In the spin-blockade regime, the double-dot is tuned 
such that one electron always resides in the right dot, 
and a second electron can tunnel from the left reservoir 
through the left and right dots, to the right reservoir [3]. 
This current-carrying cycle can be described with the 
occupations (m, n) of the left and right dots: (1, 1) — > 
(0,2) — > (0,1) — > (1,1). When an electron enters the 
left dot and forms a double-dot singlet state Sn with the 
electron in the right dot (S = | fl) — | J.|), normalization 
omitted for brevity), it is possible for the left electron 
to move to the right dot, because the right dot singlet 
state S02 is energetically accessible. Next, one electron 
tunnels from the right dot to the right lead and another 
electron can again tunnel into the left dot. If, however, 
the two electrons form a double-dot triplet state Tn, the 
left electron cannot move to the right dot, as the right 
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dot triplet state T02 is much higher in energy (due to 
the relatively large exchange splitting in a single dot). 
The electron can also not move back to the lead and 
therefore further current flow is blocked as soon as any 
of the (double-dot) triplet states is formed (Fig. la,b). 

Spin blockade only occurs if at least one of the eigen- 
states of the system Hamiltonian is a pure triplet state. 
If processes are present that induce transitions from all 
the three (1,1) triplet states to the (1,1) singlet state, 
spin blockade is lifted and a current will flow. As we will 
see below, the presence of the nuclear spins in the host 
semiconductor can give rise to these kind of transitions. 
This can be seen most easily by adding the effect of the 
hyperfme interaction to the system Hamiltonian. 



C. System Hamiltonian 

The system Hamiltonian is most conveniently writ- 
ten in the two-electron singlet-triplet basis with the 
quantization-axis in the z-direction. The basis states are 
Six, 1\X, T{1, T°! and £02 • The subscript m, n denotes 
the dot occupancy. We exclude the T 02 state from the 
model, because this state is energetically inaccessible and 
therefore does not play an important role in the trans- 
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FIG. 1: a) A schematic of the double dot and the electro- 
chemical potentials (energy relative to the (0,1) state) of the 
relevant two-electron spin states. For Alr > t, transitions 
from the S11 state to the S02 state are possible via inelastic 
relaxation with rate Ti„. Spin blockade occurs when one of 
the Tli states is occupied, b) Similar schematic for Alr — 0, 
where the singlet states are hybridized. Also in this case, 
spin blockade occurs when one of Tfj states is occupied, c) 
Energy levels as a function of detuning. At Alr = 0, the 
singlet states hybridize into bonding and anti-bonding states. 
The splitting between the triplets states corresponds to the 
Zeeman energy gnBB ext - 



port cycle. Furthermore, we neglect the thermal energy 
kT in the description, which is justified when the bias 
over the two dots is much larger than kT. The system 
Hamiltonian is given by 

H = - A Lfl |5 02 )(S*o2|+t(|5ii)<5*o2| + |5*o 2 )<5ii|) 
- g» B B ext (|T+) (T+ 1 - |T n ) (Tu \), (1) 

where Alr. is the energy difference between the \Sn) 
and IS02) state (level detuning, see Fig. la), t is the tun- 
nel coupling between the |Sn) and |5o2) states, B ext the 
external magnetic field in the z-direction and ££(#) the 
spin operator along z for the left (right) electron. The 
eigenstates of the Hamiltonian (1) for finite external field 
are shown in figure lc. For Alr < t, the tunnel coupling 
t causes an anti-crossing of the |<Sii) and |<So 2 ) states. For 
Al;? < 0, transport is blocked by Coulomb blockade (i.e. 
the final state |5o 2 ) is at a higher energy than the initial 
state S11). For Alr > 0, transport will be blocked when 
one of the three triplet states becomes occupied (spin 
blockade). In Fig. la and b, we distinguish two regimes: 
Alr > t where the (exchange) energy splitting between 
Tj^ and Su is negligibly small and transitions from Su 
to S02 occur via inelastic relaxation with rate Ti n and the 
energy. A different regime holds for Alr < t, where Su 
is coherently coupled with S02 giving rise to a finite (ex- 
change) splitting between and the hybridized singlet 
states. We will return to this distinction in the discussion 
below. 



D. Singlet-triplet mixing by the nuclear spins 

The effect of the hyperfme interaction with the nuclear 
spins can be studied [12] by adding a static (frozen) ef- 



fective nuclear field B^, 
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For the sake of convenience, we separate the inhomoge- 
neous and homogeneous contribution, for reasons which 
we will discuss later. Considering the nuclear field as 
static is justified since the tunnel rates and electron spin 
dynamics are expected to be much faster than the dy- 
namics of the nuclear system [10, 13, 14]. Therefore, 
we will treat the Hamiltonian as time- independent. The 
effect of nuclear reorientation will be included later by 
ensemble averaging. 

We will show now that triplet states mix with the Su 
state if the nuclear field is different in the two dots (in 
all three directions). This mixing will lift spin blockade, 
detectable as a finite current running through the dots for 
Ala > 0. The effective nuclear field can be decomposed 
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FIG. 2: a) Observed current flow in the inelastic transport 
regime (gfiB^LR 2> t) due to singlet-triplet mixing by the 
nuclei, b) Electrochemical potentials in the presence of H nuc i 
(t ~ ABjv). Singlet and triplet eigenstates are denoted by 
red and blue lines respectively. Hybridized states (of sin- 
glet and triplet) are denoted by dotted purple lines. For 
gfiBBext S> t, g/isABjv, the split-off triplets (T-j^ andTj~[) are 
hardly perturbed and current flow is blocked when they be- 
come occupied. Parameters: t = 0.2/^eV, g/isBjv,L = (0.1,0,- 
0.1)/xeV, gHBB NtR =(-0.1 ,-0.2,-0.2)peV and gf i B B ext =2lJ.eV. 

in a homogeneous and an inhomogeneous part (see right- 
hand side of (2)). The homogeneous part simply adds 
vectorially to the external field B ext , changing slightly 
the Zeeman splitting and preferred spin orientation of the 
triplet states. The inhomogeneous part AB^r = Bj^ — 
B^ on the other hand couples the triplet states to the 
singlet state, as can be seen readily by combining the 
spin operators in the following way 

SI- SI = ^{\ Sll ){T^\-\S lx ){T+\ +h.c.) 

SI- SI = ^(ilSnKT^l-ilSnXT+l + h.c) 

Sl~S z R = l(\S 11 ){T° 1 \ + \T^ 1 ){S 11 \). (3) 

The first two expressions reveal that the inhomogeneous 
field in the transverse plane ABff, &-B v N mixes the \T^[) 
and |Tj^) states with the \Sn). The longitudinal com- 
ponent AB^ mixes IT^) with |5n) (third expression). 
The degree of mixing between two states will depend 
strongly on the energy difference between them [5], In 
the case of g/iBB ext ,t < g/iBy/ (AB'fj), the three triplet 
states are close in energy to the |5n) state. Their 
intermixing will be strong, lifting spin blockade. For 
gfi B B ext > t,gfi B \/( AB N) the \ T ii) and \ T n) states 



are split off in energy by an amount of g[i B B ext . Con- 
sequently the perturbation of these states caused by the 
nuclei will be small. Although the iT-j^) remains mixed 
with the 15*11) state, the occupation of one of the two 
split-off triplet states can block the flow through the sys- 
tem. 

The effect of nuclear mixing is shown in Fig. 2 [5]. 
The observed current flow through the system is typi- 
cally in the order of a few hundreds of fA (Fig. 2a). 
At zero field, where the mixing is strongest, the current 
flow is largest. Increasing the field gradually restores 
spin blockade. Fig. 2b shows the energy levels for zero 
and finite external field. The theoretical calculations of 
the nuclear-spin mediated current flow (obtained from 
a master equation approach) are discussed in references 
[12, 15]. 



E. Oscillating magnetic field and rate equations 

So far, we have seen that the occurrence of transitions 
between singlet and triplet spin states are detectable as 
a small current in the spin blockade regime. We will 
now discuss how this lifting of spin blockade can also be 
used to detect single spin rotations, induced via electron 
spin resonance. The basic idea is the following. The 
basic idea is the following. If the system is blocked in 
e.g. |T) IT) i an d the driving field rotates e.g. the left 
spin, then transitions are induced to the state ||)|T)- 
This state contains a singlet component and therefore a 
probability for the electron to move to the right dot and 
right lead. Inducing single spin rotations can therefore 
lift spin blockade. 

However, together with the driving field, the spin tran- 
sitions are much more complicated due to the interplay of 
different processes: spin resonance of the two spins, inter- 
action with the nuclear fields, spin state hybridization by 
coherent dot coupling and inelastic transitions from the 
S(l,l) state to the S(0,2) state. In order to understand 
the interplay of these processes, we will first model the 
system with a time-dependent Hamiltonian and a den- 
sity matrix approach. Next, we will discuss the physical 
interpretation of the simulation results. 

The Hamiltonian now also contains a term with an os- 
cillating magnetic field in the x-direction with amplitude 
B ac 

H ac (t) = g[i B B ac sm(uT) ■ (S X L + S%). (4) 

We assume that B ac is equal in both dots, which is a 
reasonable approximation in the experiment (from simu- 
lations we find that the difference of B ac is 20% at most 
[1]). We assume B ext ^> Biq,B ac , which allows applica- 
tion of the rotating wave approximation [16]. Therefore, 
we will define B\ = ^B ac , which is in the rotating frame 
the relevant driving field for the ESR process. 

In order to study the effect of ESR and the nuclear 
fields that are involved in the transport cycle, we will 
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construct rate equations that include the unitary evo- 
lution of the spins in the dots governed by the time- 
dependent Hamiltonian. This approach is based on the 
model of reference [12], where the Hamiltonian contained 
only time-independent terms. Seven states are involved 
in the transport cycle, namely the three (1,1) triplets 
I7W), the double and single dot singlet states |Sn) and 
\So2) and the two (0,1) states | foi) and I loi), making 
the density operator a 7 x 7 matrix. The rate equations 
based on the time-independent Hamiltonian are given in 
[12]. These are constructed from the term that gives the 
unitary evolution of the system governed by the Hamilto- 
nian (H = H a + H ac ) dp k /d,T = -^{k\[H, p]\k), together 
with terms that account for incoherent tunnelling pro- 
cesses between the states. The rate equations for the 
diagonal elements are given by 
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The rate equations for the off-diagonal elements are 
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FIG. 3: Time evolution of the diagonal elements of the density 
matrix for one particular nuclear configuration. Parameters: 



100 mT, B 



SMslOOmT, B ext 

Si = 1.3 mT 



B§, x , y , z =(0,0,0) 
MHz, hV in = gp B B^ z and A L r 



N,x,y,z =(0,0,2.2) mT, 
T L = 73 MHz, T R = 73 
200/ieV, i=0.3 ^eV. 



dpjk 
dr 



{j\[H,p}\k)--{^+T k )p jk 



(6) 



where the indices j, k € {T 11; Su, S02, Toij I01) label the 
states available to the system. The tunneling/projection 
rates r j equal T in and for the 1511) and IS'02) states 
respectively, and equal zero for the other 5 states. The 
first term on the right-hand side describes the unitary 
evolution of the system, while the second term describes 
a loss of coherence due to the finite lifetime of the sin- 
glet states. This is the first source of decoherence in our 
model. The second one is the inhomogeneous broaden- 
ing due to the interaction with the nuclear system. We 
do not consider other sources of decoherence, as they are 
expected to occur on much larger timescales. 

Because we added a time-dependent term to the 
Hamiltonian (the oscillating field), we numerically calcu- 
late the time evolution of p(t), treating the Hamiltonian 
as stationary on the timescale Ar <^2it/uj. To reduce the 
simulation time, we use the steady state solution p T ^oo 
in the absence of the oscillating magnetic field as the ini- 
tial state p(r — 0) for the time evolution. At t = the 
oscillating field is turned on and the system evolves to- 
wards a dynamic equilibrium on a timescale set by the 
inverse of the slowest tunnelling rate T. This new equi- 
librium distribution of populations is used to calculate 
the current flow, which is proportional to the occupation 
of the IS'02) state (/ = eTnps 02 ). An example of the time 
evolution of the density matrix elements is shown in Fig. 
3. The figure clearly reveals that the blockade is lifted 
when the oscillating field is applied. This is visible as an 
increase of the occupation of the IS'02) state. 

In order to simulate the measured current flow we have 
to consider the fact that the measurements are taken 
with a sampling rate of 1 Hz. As the timescale of the 
nuclear dynamics is believed to be much faster than 1 
Hz [10, 13, 14], we expect each datapoint to be an in- 
tegration of the response over many configurations of 
the nuclei. The effect of the evolving nuclear system is 
included in the calculations by averaging the different 
values of the (calculated) current flow obtained for each 
frozen configuration. These configurations are randomly 
sampled from a gaussian distribution of nuclear fields in 
the left and right dot (similar as in [12]). Because the 
electron in the two dots interact with different nuclear 
spins, the isotropic gaussian distributions in the two dots 
are uncorrelated, such that yj (AB^) = \f2\J (B 2 N ) and 
(B%,x) — (B%,y) ~ (B%,z)- F° r the sake of convenience 
we define 
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Simulation: Inelastic transport (A LR »t) 
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Simulation: Resonant transport (A LR =0) 
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FIG. 4: a). Calculated average current flow in the inelastic 
transport regime. Parameters: fioj = g/iBlOOmT, B ea; t =100 
mT, ctjv =2.2mT, Si = 1.3mT, F L m = 73 MHz, i=0.3 /xeV 
and A_lh = 200 /^eV. Results are similar for any value for t, 
provided that Alb. 3> t. b) Calculated average current flow 
in the resonant transport regime at zero detuning for differ- 
ent values of t. Parameters: hu = g/iBlOOmT, on =2.2mT, 
Bi = 1.3mT, F L , R = 73 MHz, F ln = and A LR = 0. Av- 
eraged over 400 nuclear configurations for t/(gfiBO~N) > 0.5 
and 60 configurations for t/(gfiB&N) =0.5. Simulation car- 
ried out for positive magnetic fields only; values shown for 
negative fields are equal to results obtained for positive field, 
c) Experimental data from Ref [1] with (curve offset by 100 
fA for clarity) and without oscillating magnetic field. The fre- 
quency of the oscillating magnetic field is 460 MHz and the 
applied power is -16dBm. Simulation carried out for posi- 
tive magnetic fields only; values shown for negative fields are 
equal to results obtained for positive field. 



F. Simulation results and physical picture 

An example of the calculated (average) current flow 
as a function of B ext (Fig. 4a, b) shows a (split) peak 
around zero magnetic field and two satellite peaks for 
B ext = ±huj/(gfj,B), where the spin resonance condition 
is satisfied. This (split) peak at B ext — is due to singlet- 
triplet mixing by the inhomogeneous nuclear field, and 
the splitting depends on the tunnel coupling, similar as 
the observations in [5]. The response from the induced 
spin flips via the driving field is visible for the both in- 
elastic and resonant transport regime, and the current 
flow has comparable magnitude to the peak at B ext — 0. 
The satellite peaks are also visible in the experimental 
data from [1] (also shown here in Fig. 4), although the 
shape and width of the satellite peaks arc different, as 
we will discuss later. 

We want to stress that the ESR satellite peaks only 
appear when an inhomogeneous nuclear field is present 
in the simulations. In other words, for ABn — and 
B\ equal in both dots, spin rotations are induced in both 
dots at the same time and at the same rate. Starting, 
for example, from the state \T^[) = | | j) transitions are 
induced to the state | J, J,) via the intermediate state | | 
+ DIT + |)/\/2= (|T+) + |T 1 - 1 )+2|T 1 ° 1 ))/x/2. No mixing 
with the singlet state takes place (the evolution is in the 
triplet-subspace) and no current will therefore flow. 

The ESR sattelite peaks are visible for both resonant 
and inelastic transport regime (Figs. 4a, b). For the res- 
onant transport regime, we see that for tja-^ < 5 the 
sattelite peak increases in height when increasing t, sim- 
ply because the coupling between the two singlet states 
increases. However, further increasing t reduces the sig- 
nal, and this is because the exchange splitting then plays 
a more important role. Namely, increasing the exchange 
splitting reduces the mixing between the state with 
the hybridized singlet state by the nuclear field gradi- 
ent. This mixing is a crucial element for detecting the 
induced rotations of one of the two electron spins. In 
the inelastic transport regime, this exchange splitting is 
negligibly small and therefore the height of the sattelite 
peak depends only on Ti n and the driving field B±. 

A study of the height of the satellite peak as a function 
of B\ reveals a non-monotonous behaviour, which can 
be seen in Fig. 5a. The physical picture behind this 
behavior is most easily sketched by distinguishing three 
regimes: 

1. For B\ < <Jn,zi for most of the nuclear configu- 
rations the spin in at most one of the two dots 
is on resonance, so spins are flipped in cither the 
left or right dot. In that case transitions are in- 
duced from e.g. | |t) to | Tl> = Pn) + fT^) or 
I IT) = Pii) — Pn). The resulting current flow ini- 
tially increases quadratically with B\ , as one would 
normally expect (Fig. 5a). 

2. For B\ 3> (Jn.z, for most of the nuclear configu- 
rations two spins are rotated simultaneously due 
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FIG. 5: Height and width of the ESR satellite peak, a) Cir- 
cles: calculated ESR peak height as a function of driving am- 
plitude B\. Parameters: hui = g^islOOmT, B ex t =100 mT, 
crjv =2.2mT, T L>R = 73 MHz, t=0.3 fieV, hV in = gfi B a N 
and Alr = 200 fieV. Lines are the current measurements 
for 2 different values of Alr- The measurements show time- 
dependent (telegraph type) behavior. Therefore, the curves 
are obtained by repeating sweeps of B± and then selecting 
the largest current value for each value of B\. b) Calculated 
width of the ESR satellite peaks as a function of B\. For 
small ESR power the peak is broadened by the random nu- 
clear fluctuations, at high powers it is broadened by B\. 



to power broadening of the Rabi resonance. The 
stronger B\ , the more the transitions occur only in 
the triplet subspace (the driving field B\ that ro- 
tates two spins dominates the S — To mixing by the 
nuclear spins). As a result, the current decreases 
for increasing B\. 

3. If B\ ~ <jn,z the situation is more complex because 
both processes (rotation of 2 spins simultaneously 
and transitions from to 6*11 ) are effective. We 
find that if both processes occur with comparable 
rates, the overall transition rate to the singlet state 
is highest. This is the reason why the current has 
a maximum at B\ m <jn,z (Fig- 5a). 

The experimental data of the ESR satellite peak height 
(normalized by the zero-field current flow) for two dif- 
ferent values of A l r are shown in Fig. 5a. In order to 
compare the experimental results with the model we have 
estimated the rate from the measured current flow at 
B e xt = (we found similar values for both curves). The 
agreement of the experimental data with the model is 
reasonable, as it shows the expected quadratic increase 
with j5i, as well as a comparable peak height. However, 
we see that variations of the level detuning A^a can re- 
sult in considerable differences of the measured ESR peak 
height. We have two possible explanations for the devi- 



ations of the experimental data with the model. First, 
we have found experimental signatures of dynamic nu- 
clear polarization when the ESR resonance condition was 
fulfilled. We expect that this is due to feedback of the 
electron transport on the nuclear spins (similar to that 
discussed in [11, 15, 17]), although the exact processes 
are not (yet) fully understood. Second, unwanted electric 
fields affect the electron tunnelling processes, but are not 
taken into account in the model. We expect that these 
electric fields will not change the location and width of 
the ESR sattelite peaks because this field does not couple 
the spin states. It is however possible that the height of 
the satellite peak is altered by the electric field because if 
can affect the coupling between the S(0,2) with the S(l,l) 
state. 

Finally, we discuss the width of the ESR satellite peak 
(Fig. 5b). If the inelastic tunnelling process between 
the dots (with rate Y^ n ) and B\ are both smaller than 
cjv.z, the ESR peak (obtained from simulations) is broad- 
ened by the statistical fluctuations of the effective nuclear 
field. For high Bi, the width approaches asymptotically 
the line with slope 1 (see Fig. 5b). In this regime, the 
peak is broadened by the RF amplitude B\. In the ex- 
periment [1] , the shape of the satellite peak was different 
(flat on top with sharp edges) than expected from the 
model. Furthermore, the FWHM was larger than ex- 
pected from just (Jn,z- We attribute this to feedback of 
the ESR-induced current flow on the nuclear spin bath. 
As a result, a clear FWHM increase with B\ could not 
be observed. 

It should be noted that in the simulation the central 
peak is broader than the satellite peaks. From study- 
ing the influence of various parameters in the model, we 
conclude that the greater width of the central peak is 
caused by the tranverse nuclear field fluctuations (Bjq,x 
and Bff, y ), which broaden the central peak but not the 
ESR satellite peaks. 

We conclude that the model discussed here qualita- 
tively agrees with the main features that were observed 
in the double dot transport measurements that aims at 
detecting (continuous wave) ESR of a single electron spin. 
The details of the ESR satellite peak height and width 
do not agree quantitatively with the model. We believe 
these deviations can be attributed to unwanted electric 
fields and feedback of the electron transport on the nu- 
clear spin polarization. Improving the understanding of 
these feedback mechanisms remains interesting for future 
investigation as it might point towards a direction to mit- 
igate the decoherence of the electron spin [12, 18]. 
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